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Abstract 



We consider a system of spinless fermions in a nearly one-dimensional cylindrical trap. 
We introduce four-fermion interaction terms through which particles in low-lying modes 
of the trap may interact with particles in radially or angularly excited states. Treating 
these interactions as perturbations about a purely one-dimensional Luttinger liquid, we 
calculate the zero-temperature energy shifts of the bosonic quanta of the Luttinger model. 
We perform this calculation is two ways. First we use the bosonized forms of the fermion 
operators. Then we introduce a simplified, physically-motivated model; in this model, the 
methods used immediately generalize, so that we may also calculate the energy shifts of 
the Tomonaga-Luttinger model. By equating the results from the two disparate models, 
we obtain the numerical value of the Luttinger liquid bandwidth for this system. 
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1 Introduction: fermions trapped in one dimension 



A low-temperature system of fermions in one dimension has the known property of 
bosonization PH El IS] ■ The coherent fermion-hole excitations of the system behave like 
bosons, with a spectrum that remains exactly solvable even in the presence of certain arbi- 
trarily strong interactions. This phenomenon is described by the Luttinger and Tomonaga- 
Luttinger models. 

The main physical applications of the Luttinger model occur in the context of the quan- 
tum Hall effect [3J |5] . However, there is another situation in which we expect Luttinger- 
type behavior to arise quite naturally. If a large number of degenerate fermions are 
trapped in a potential well with a very high aspect ratio, so that the particles are effec- 
tively constrained to move along a single axis, we would expect one-dimensional behavior 
to dominate. We shall consider a system of this type. Since the trapping potential cannot 
be purely one-dimensional, the second and third dimensions will affect the energy levels 
of the system. We shall calculate these energy shifts, by treating the effects of the small 
transverse dimensions as perturbations around a purely one- dimensional model. Since the 
field of fermion trapping is advancing very rapidly |H1 13 El E] , such calculations should be 
highly relevant in the analysis of future research. 

After a brief discussion of trapped fermions and the one-dimensional Luttinger model, 
we shall introduce a physically motivated interaction Hamiltonian in Section El We then 
calculate the perturbative energy shift due to this interaction in two different ways. In 
Section 01 we use bosonized forms for the fermion operators. The results of the bosonized 
calculation depend upon the cutoff of the model. We may determine the physically mean- 
ingful value of the cutoff by comparing the bosonized result with the result of our second 
calculation (Section 0J), which uses an entirely different, simplified model of the one- 
dimensional Fermi gas. Finally, in Section we extend our results to systems with 
additional interactions and discuss the significance of our findings. Our main result is 
expression ([55)1. which describes the energy shifts of the bosonic modes. 

We shall consider a system with a large number of fermions trapped in a cylindrical 
potential. In the cylidrical coordinates, (z,8,p), the wave function is periodic in z with 
period L and has a Dirichlet boundary condition at p = b <C L. The boundary condition 
for p corresponds to a trap with hard walls, but our general treatment will also be appli- 
cable to other forms for the trapping potential. The eigenfuctions of this trap are (setting 
h = l) 

e ikz e m Jt(ji,n P p/b) 

Uk £ n p — — 7= — 1= 1 i (1) 

VLV2nbJ-Jt +1 (jt, np )Je-i{jt,n P )/2 

where k = ^ for some integer r, £ is an integer, and n p is a positive integer, which 
represents the number of radial nodes in the wave function. The number ji tJl is the n p -th 
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zero of the Bessel function Jg. The energies of these modes are 



Jl,n p 



(2) 



so there are energy gaps A^ np = ^ t ' np \ m y2°' " > betwen the n p = £ = states and the higher 
exicited states. 

In this paper, we are interested in the situation where the Fermi level is low enough 
so that only the n p = £ = states are occupied, and the exicted states of the trap can 
enter only virtually. This requires that the transverse size b of the trap be of the order 
of the mean spacing between adjacent particles or smaller, (and most of our calculations 
will also require that the temperature of the system be vanishing). When no radially or 
angularly excited states are occupied, k becomes the only relevant quantum number, and 
the system becomes effectively one-dimensional. This reduction in dimensionality allows 
us to make use of the peculiar theory of one- dimensional degenerate fermions, with its 
property of bosonization. 

The phenomenon of bosonization is expressed in the Luttinger model, in which the 
noninteracting (spinless), one-dimensional fermionic Hamiltonian is approximated by 

H = v F J2 [( k ~ k F )a\a k + {-k - k F )b\b k , (3) 

k 

where a (b) operators correspond to right- (left-)moving fermions, and v F is the Fermi 
velocity. The linearization of the energy spectrum comes from the fact that the energy 

required to excite a particle near the Fermi surface is ^t,^ = Vp( l + 2m> ^ the 

excitation momentum q is small compared to the Fermi momentum, them we may neglect 
the q 2 term. In the Luttinger model, the further approximation is made that this Hamil- 
tonian is valid for all momenta k, but that states corresponding to a k and b- k for k < 
are filled [T0| lllj. The introduction of fictional negative-energy states, which lie far from 
the Fermi surface, should have little effect on the physics at low temperatures. 
With these approximations, the Hamiltonian may be expressed in the form 

h* = ^EM!)p+N+/t(-W?l (4) 



g>0 



where the p± operators 



p+(q) = Y. a l+q a k ( 5 ) 

k 

p-(q) = EO* ( 6 ) 
k 

obey the Bose commutation relations [p + (—q), p + (q')] = [p-(q), p~(— q')] = &qq'^ and 
[p + (q), p~(q')} = 0. The bosonic quanta created and annihilated by the p± operators have 
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energies uj q — Vp\q\. It is also possible to express the fermion field operators in terms of 
the bosonic p± operators. We shall make use of this representation in Section El 

The Tomonaga-Luttinger model is formed by adding scattering effects [TJ], through 
the interaction Hamiltonian 

R l = ^E{ 2 ^(g)p+(g)P-(-?) + 9A{q)[p + {q)p + {-q) + P-(q)p-(-q)]}. (7) 

The #2 and g& are the Fourier transforms of the interaction potentials. This interaction 
(like the similar interactions to be introduced in Section |2J) describes forward scattering 
processes between fermions, with momentum exchange q. The Tomonaga-Luttinger model 
is exactly solvable, since H + Hi remains bilinear in the p±. We shall return to this 
interaction in Section but first we shall consider interactions between the Luttinger 
liquid bosons and fermions in excited states of the trap. 



2 Interaction Hamiltonians 

The Luttinger model should be a good description of a system of fermions in a narrow, 
axially symmetric trap, in which no radially or angularly excited states are occupied. We 
shall add to this system additional two-fermion interaction terms. These new interactions 
will be similar in structure to Hi (when Hi is expanded in terms of fermion operators). 
These interactions must conserve linear momentum along and angular momentum around 
the axis of the trap, as well as fermion number. One such interaction term is 

H' = ^2 C k ,k>,qC k - q dk{a k , +q ak> + b k , +q b k/ ). (8) 

k,k',q 

The operators c k and d k correspond to arbitrary fermion states, subject to the condition 
that at least one of them corresponds to an excited state of the trap (i.e. is not an a k or 
a b k ). This ensures that we shall be studying new effects. The effective coupling C k) k',q 
will depend upon the radial and angular quantum numbers corresponding to c k and d k 
(in part because C kj k',q includes integrals over the radial and angular eigenmodes of the 
trap). Angular momentum conservation demands that c k and d k annihilate states with 
the same angular momentum. 

We are only interested in interactions that will renormalize the bosonic energies. The 
simplest way for bosonized operators to arise from from an interaction Hamiltonian such 
as H 1 is for C ktk r >q to be effectively independent of k and k' . Then the a k , +q and a k > 
operators combine to form a boson operator, as in (JHJ), and the b k , + and b k i combine 
similarly. This leads to an interaction of the form 

H * = ^5>(<?) (z4- 9 4) \p + (q) +P-{q)] +h.c. (9) 

q^O \ k I 
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We shall calculate the energy shifts of the low-lying bosonic states to second order in 
H 2 . (The first-order shifts obviously vanish.) If both c k _ q and d k correspond to excited 
states of the trap, then the operator H 2 preserves the number of excited fermions. Hence, 
H 2 has zero matrix element between low-lying (Luttinger liquid) states and states with 
excited fermions. So we only get nonvanishing contributions from terms where one of c k 
and d k corresponds to an excited fermion state, while the other corresponds to a low-lying 
mode. Moreover, the excited mode must have angular momentum t = about the axis 
of the trap. 

Let us fix q > 0. We shall calculate the energy shift for the p+(q) quanta. In order for 
these quanta to be well-behaved Luttinger liquid phonons, we must have u q <C ep = ^k F , 
or q <C hp. We may write the interactions affecting the p+(q) quanta as 



H. 



2, +q 




E 

. k 



(4-<A + d\ c k ) 



p+(q) 



. k 



(cl+Jk + d\ + c k ) 



P+(-q)\, (10) 



where d k annihilates an excited fermion and c k annihilates a low-lying fermion. [If c k 
annihilates a right-moving fermion, additional terms are actually be required in order 
make H 2+q hermitian. However, these terms do not involve p + , so their matrix elements 
are smaller by a factor of 0(L~ l l 2 ) and may be neglected.] We shall denote by A the 
energy gap A 0jn „ for the particular £ = states corresponding to the operators d k . 



3 Energy shift calculated with bosonized operators 

We shall calculate the perturbative energy shift arising from (fTUj) in two different models. 
In this section, we use the bosonized field operators of the Luttinger model to perform the 
calculation. It is possible to obtain a closed-form expression for the energy shift in two 
different situations — when the temperature vanishes or when the energy gap is very large. 
(In each of these cases, the "transverse temperature," which is the effective temperature 
for excitations of the trap's radial and angular modes, is zero.) In these two situations, 
we find that the leading contributions to the the infinite sum over intermediate states 
may be resummed. This resummation allows us to find a simple expression for the energy 
shift. 

We shall restrict our attention to the case in which c k is a k (rather than b k ). In this 
case, the matrix element which we must calculate is ({m +p }; K\H 2 , +q \{n +p }). The state 
is a number state of the Luttinger model; the {n +p } are the occupation numbers 
of the p + quanta. The state |{m +p }; K) has differing occupation numbers {m +p }, as well 
as an excited state fermion with momentum K. 

The matrix element is simply 

({m +p }; K\H 2)q \{n +p }) = -g(q)({m +p }] K\ £ [d\_ q a k p + {q) + d\ +q a kP+ (-q)\ \{n +p }) 
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L k ' V27T 



+29(Q)J2 5 K,k+q]J (j^ n +q {{m +p }\a k \{n'[ p }), (11) 

where the occupation numbers {n' +p } and are the same as the {n +p }, except that 

n' +q = n +q + 1 and = n +q — 1. As L — ► oo, {ra+ p } and {ri" p } differ from {n +p } 
only on a set of measure zero, and we may replace the matrix elements ({m +p }\ak\{n' +p }) 
and ({m +p }\a k \{n'[ p }) appearing in (|TTJ) by ({m^}!^!^^}) and ({m'^j-la^fn+p}}, 
where m' +(J = m +q — 1 and = m +q + 1. We shall eventually be summing over all 
possible sets of occupation numbers {m +p }, so we may ultimately replace both matrix 
elements in (JTTJ) with {{m +p }\ak\{n +p }) . So we see that the crucial quantity to calculate 
is ({m +p }\ak\{n +p }) , for an arbitrary set of occupation numbers {m +p }. We shall perform 
this calculation by taking a Fourier transform, 

1 r L 

({m +p }\a k \{n +p }) = -= / dx e~ lkx ({m +p }\^ + (x)\{n +p }) , (12) 



VZJo 

where ip + is the field operator for the right-moving fermions (while tp_ is the left-moving 
fermions' field operator). So we must now turn our attention to the operators ip±. 
The bosonized form of the fermion field operator is ^3 El 113 EH fl"T| 

, 4>±{ x ) — i - U± exp[±ik F x =F i(j)[x) + iQ(x)\. (13) 
V27ra 

In the pure Luttinger model, this is only an operator identity when we take the limit 
a — > [THj. However, we shall interpret a as a finite cutoff, so that vpa' 1 is the bandwidth 
of the Luttinger liquid [2J EU EH I2H HI] • The finiteness of a corresponds to a deviation of 
the energy spectrum from the Luttinger form far from the Fermi surface. We shall discuss 
which values of a are physically meaningful in Sections |U and El The U± operators decrease 
the total particle number on their corresponding branches by unity; these operators are 
necessary because the bosonic fields conserve the total fermion number. 

The bosonic quantities <fi (the boson field) and 9 (the integral of the momentum con- 
jugate to 0) are given by 

-{(/, = -t E -e- alpl/2 e- ipx [p + (p) + p_(p)] + m N+ + N ~ x (14) 
L p%P ' L 

16 = -e- a ^ 2 e- ipx \p + (p) - p(p)] +tir N+ ~ N ~ Xj (15) 

L p7oP ' ' L 

where N + and iV_ are the excess numbers of right- and left-moving fermions relative to 
the ground state, respectively. So the argument of the exponential in ip± is 

± ikpx T i<p + i9 = ±ik F x T ^ E -e- alpl/2 e- ipx p±(p) ± ^^x. (16) 
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For states sufficiently close to the ground state, we may neglect the N± terms; this will 
be equivalent to negelcting an term in the argument of a momentum-conserving 

^-function. So we have 

± ikpx T i(j) + i9 = ±ik F x T^T -e~ ap/2 \e~ ipx p± (p) - e ipx p±(-p)} . (17) 

L ^>P 1 J 

Since [p±(p), p±{p')] — if p ^ ±p', the elements of the sum in (fP7|) all commute with 
one-another. This allows us to write ip± in the product form 

^ ± (x) = - 1 L=U±e ±lkFX TT exp j± [^] e~ ap/2 \e ipx p±(-p) - e- ipx p±(p)} 1. (18) 
V2^o7 { \pLj L Jj 

So the study of the matrix elements of ip± between different Luttinger liquid states leads 
naturally to the study of harmonic oscillator matrix elements of the form (m|e AAt ~ A *" 4 |n), 



where and A are raising and lowering operators. In our case, since ^2ir/(pL)p± is a 
canonically normalized ladder operator, 

A = T J^re Tipx e- ap/2 . (19) 



It is important to note that A is 0(L 1 ^ 2 ). 

In terms of / = m — n, the general formula for the matrix element in question is 

1 /ml\^ sgn(0 

(m\e XAt - x " A \n) = ± IHt) A CW+0/a(_ A *)(l«l-0/a 
\l ! \ n\ I 



{e |A|2/2 F [max(m, n) + 1; |/| + 1; -|A| 2 ]} , (20) 



where F(a; b; z) is a confluent hypergeometric function. This formula is derived in Ap- 
pendix 1X1 

The formulas above pertain to both ip + and ip-, but we shall henceforth restrict our 
attention to ip + only. We shall now proceed to calculate the energy shift for the case in 
which the temperature is effectively zero when compared to the transverse energy scale 
of the trap. The fact that the temperature vanishes allows us to make some substantial 
simplifications of the mathematics. We shall show explicitly that these simplifications are 
not valid (and would give rise to unphysical results) when the "transverse temperature" 
is nonzero. 

It may initially appear that the bracketed factors in (|20j) may be neglected, since both 
e' A ' I 2 and F[max(m, n) + 1; |/| + 1; — |A| 2 ] are power series in |A| 2 , and |A| 2 is C(L _1 ). 
Each power series is dominated by its first term in the large L limit, and for both series, 
that term is unity. However, these terms may not actually be neglected in the L — » oo 
limit, because there are an infinite number of such terms, corresponding to all the allowed 
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momenta. Since the density of states in momentum space is proportional to L, the infinite 
product does not approach unity. 

Let us denote the infinite product in question by V; that is, V is defined to be 



V 



II exp 

p>0 









F 







max (m +P7 n 




(21) 



When we calculate the energy shift, we sum up terms consisting of a matrix element 
squared divided by an energy difference. If the characteristic energy differences are O(L), 
then the sum of the matrix elements squared must also be O(L) if we are to obtain a 
nonvanishing contribution. This is only possible if the occupation numbers of a macro- 
scopic number of states states are O(L). This is a natural situation at finite temperature; 
however, at T = 0, the occupation numbers must be small, so that the excitation energy 
of the system is not extensive. 

So at T = we only get a nonvanishing contribution when the difference in energies 
between \{m +p }) and grows more slowly than L. That is, almost all of the m +p 

must be equal to the n +p . We may then separate out those values of p at which m +p and 
n +p differ and get 



V = J] exp 

x n \ ex p 

p,l +p =0 I 




max (m +p , n +p ) + 1; \l +p \ + 1; - (^j^J e ap 



n +p + !;!;- fej e 



-ap 



(22) 



The number of terms in the first product must grow more slowly than L. Since each of 
these terms has the form 1 + 0(L _1 ) as L —>■ oo, the entire product is one in that limit. 

The second product in (}2*2*|) is more complicated. Using the Kummer transformation 
formula F(a; b; z) = e z F{b — a; b; —z) and the fact that F(—n; 1; z) is just the Laguerre 
polynomial L n (z) (221 > we ma y rewrite V as 



V= I] exp 

p,l +p =0 I 




J n +V 




-ap 



(23) 



We may evaluate this explicitly at T = 0, where n +p = for almost all p. The product 
of the Laguerre polynomials is one, and all that remain are the exponentials. Each 
momentum p represented in the product has the form p(r) = 2-Kr/L, for some positive 
integer r, so we have 



V{T = 0) = exp 



E 



r J+p(r)— 



^—2-KCtr/L 

2r 



(24) 
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Since l +p = for almost all p, we may extend the sum to include all positive integers r. 
Then the sum is just the Taylor series for |log — e~ 27rQ//L ), so, to leading order in L -1 , 



we have 



P(T = 0) = *fe (25) 



The factor of y ^ will be necessary in order to make our final result finite. 

If we were to assume that the only terms which contributed had m +p = n +p for almost 
all p even when T > 0, we would obtain the expression 



v=p^\ n l 

\p,l+P=0 




(26) 



We shall use this expression to demonstrate explicitly that our assumption is inconsistent 
unless T vanishes. 

We may now turn our attention to remaining factors in the matrix elements — those 
arising from the unbracketed terms in (|20p. We shall divide the matrix elements of this 
operator into groups. Each matrix element ({m +p \\ip + \{n +p }) is labeled by an integer 
j > 0, where j is the number of quanta by which {m +p } and {n +p } differ. That is, 

j = J2\m +p -n +p \. (27) 

p>0 

According to our earlier argument, at T = 0, the relevant values of j are smaller than 
O(L), and we shall for the moment assume that this is also the case at finite T (although 
we shall eventually show that this assumption is inconsistent). 

Each of the j quanta by which {m +p } and {n +p } differ corresponds to a momentum 
Pi > 0, i G 1, 2, . . . ,j. Since j is smaller than O(L), in the continuum limit all the Pi are 
distinct, except on a set of measure zero. So we shall only need to use (j2T)j) in the case 
where |Z| < 1. Keeping this in mind, we may write the matrix elements of ip + as 

-p i f _ / l2L e -m x e - a Pi/ 2 / n _i_ i 

({m +P }\Mmn +P }) = -l==e tkFX U /fe_- aDi / 2 J- ■ (28) 



t=\ 1 J^e^e-^^p-. 



PiL 

The upper term applies if m +Pi > n +Pi (i.e. if m +Pi = n +Pi + 1), and the lower term 
otherwise. We note that this matrix element is 0{L^^ 2 ). Such a matrix element can 
only give rise to a nonvanishing energy shift if it is multiplied by an appropriate positive 
power of L. 

We may find the corresponding matrix element of by taking the Fourier transform 
(H2J- We have 



({m +p }\a k \{n +p }) = — — -= / dx exp —ikx + ik F x + V] =fip, 

V lira \ L J v L Jo \ 



x 



S 



^(-HmM'v^ 1 ■ <29) 

For each i, we have two possibilities, depending upon the relative magnitudes of m +Pi and 
n +Pi . In each case, the upper (lower) sign in the exponential corresponds to the upper 
(lower) square root. The integral over x just gives L5 k fcj?+ j-^ Tp .- The matrix element 
squared is then 

(30) 

To find the perturbative energy shift due to terms of this form, we must divide the 
matrix element squared by the energy defect and integrate over all possible momenta. 
The phase-space factors arising from the integrations will cancel out the j factors oi2n/L 
in (HUH). The energy defect is given by 



E w+p} . kTq - E {n+p} = A + — 



(kTq?-k 2 F \ ±u q + J2±u Pi - (31) 

i=l 

A + =F ?) 2 is the energy of the excited fermion; -^k 2 F is the energy lost by the 
Luttinger liquid ground state with the removal of one fermion; and ±u q + Yh=i is 
the energy change due to the phonons gained or lost in the intermediate state. There are 
2 J possible choices of signs in the sum over the quanta created or annihilated by a^. 

The momentum range over which we must integrate extends from e~ 7 p m i n = e -7 ^ 
to +00. The Euler-Mascheroni constant 7 enters when we change from a sum to an 
integral, through the formula lim m ^ +00 (j2T=i \ ~ Ji" 7") = 7- We may then absorb 
this 7-dependence into the limits of integration. When we integrate all the pi over this 
range, we must also divide by j!. This factor arises from our overcounting of intermediate 
states — permuting the Pi does not give rise to a new state. 

So the energy shift consists of terms of the form 

AEi = -ZL_i_ Wa) |» E E 1 J. (Jl)' (32) 




L6 k,k F+ j:i lTPl { n +i + 



s+l 
2 



A + h a k - s( i) 2 - k V\ + SU q + H=l ±W; 



V, 



The sum over t runs over the 2- 7 choices of sign (and of n +Pi + 1 or n +Pi ). The integrand in 
(jS!?)) behaves as p^ 1 for small values of Pi and as p^ 3 for large values of Pi, so the integral 
is dominated by the region where all the Pi are small. 

Near pi = p m i n , Pi and u p . are each O^L" 1 ), so we may neglect J2i=i TPi and Y^l=i 
unless j is O(L) or greater. Since we have assumed that j is smaller than O(L), we may 
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drop Yh=i TPi and Yl{=i ±k> Pi - We shall not drop the exp (—aY^iPi) term; this factor 
is needed to provide a cutoff for the momentum integral. 

By dropping the two sums, we simplify the integral substantially. The expression 
becomes 



V 2 



V 



EE 

s=±lt=l 



n +Pi + 1 



n 



A + h K k F ~ sq) 2 ~ k%\ + su q 



(33) 



We may now perform the sum over t, to get 
V 2 



47r 2 a 



E 

s=±l 



g(q)\ 2 q (n +q + **) 1 



oo 

Pmin/eT 



dpi 



Pi 



(2n +K + 1) 



(34) 



Furthermore, we may easily sum up all the terms with different values of j, getting the 
exponential function 



AE 4 



V 2 



s+l 



4n 2 a 



\g(q)\ 2 q (n+ q + 2 



exp 



s=±l 



(35) 



(The subscript "+" in AE + indicates that this energy shift only includes the terms arising 
from ip + — i.e. the terms in which Ck is a*..) It now remains to evaluate the integral and 
add the contributions from the terms for which Ck is bk- 

It is impossible to evaluate the integral for a general set of n +Pi . However, we may 
simplify our expression for the energy shift substantially by separating the two terms in 
the integral and transforming the integral back into an infinite sum. We have 



exp 



OO g 

dpi 



-apt 



Pi 



(2n +Pl + 1) 



exp 



e- api n+Pi 



Pi>0 



exp 



Pl >o \Pt L J 



■apt 



(36) 

The second factor on the right-hand side of (|3fij) contains the same sum that we encoun- 
tered when we evaluated V(T = 0); in fact, this term is exactly [V(T = 0)]~ 2 . So we 
define a new infinite product T> as 



V 



V exp 



n i ex p 

Pi>0 I 



+oo 



dpi 



-api 



Pi 



-ap 



'n 



+Pi 



l +Pi 



,Pi L , 



-ap. 



(37) 



T> is an overall factor multiplying our expression for AE + . It is clear that T> — 1 when 
T = (since almost all the n +p are then zero). For T > 0, T> contains all the energy 
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shift's temperature dependence, so it is natural to calculate its thermal average. We shall 
find that this average diverges for all finite temperatures; this unphysical result signals 
that we can not assume that j is smaller than O(L) unless the temperature is strictly 
vanishing. 

Denoting the inverse temperature by (3, we find that the thermal average of T> is 



{V)=M J2 {exp 

{ n + P } 



-f3v F J2 (P n + P ) 

p>0 



V 



(3* 



The sum over the {n +p } runs over all possible sets of occupation numbers, and M is the 
necessary normalization factor 



N = { Yl ex P 

A n +p} 



-f3v F J2 (p u +p) 

p>0 



n(i-e 

p>0 



-/3v F p 



(39) 



(T>) may be factorized into an infinite product of sums indexed by p. This allows us to 
calculate the thermal average separately for each value of the momentum, so (V) becomes 



(V) = TV Yl ( ex P [-P v FPn+p\ \ exp 

p>0 n+ p =0 



LVpl) e ~ aPn+p . 



n +p 



,p l , 



-ap 



■ (40) 



The sum is now in the form of a product of generating functions for (L n ) 2 . That is, 



(V) = X UJ2[L n (x)] 



2 „n 



(41) 



p>0 n=0 



with z = exp {^i) e ap ~ (3vfP 

The sum of this generating function is known [23]. In fact, when \z\ < 1, 



and x = (||J e~ ap . 



2 n 



71=0 



exp 



2xz 



2x- 



(42) 



where Jo is the modified Bessel function Io(u) = Jo(iu). The series diverges if \z\ > 1. So 
(V) converges only if exp e~ ap — PvfP^ < 1 for all p. Since (jjij e~ ap — (3vfP is a 
strictly decreasing function of p, we need only consider p = p m in- At this value of p, we 
may neglect e~ ap , so the condition for convergence becomes (3 > {ypPmin} 1 ■ As L — > oo, 
this can not be satisfied at any finite temperature. So the thermal average of T> is equal to 
unity at T = (when z = for all p) and diverges for T > 0. {For finite L and very large 
Vp, the condition that (3 > (fFPmin) -1 gives us a qunatitative estimate of how small the 
termperature must be in order for our approximation to be valid. Current experiments 
with degenerate fermions 001312] typically achieve temperatures for which f3~ l ~ e F - 
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In order for the T = approximation to be valid, the temperature must be smaller by a 
factor of O(N), so this regime is far beyond the limits of current technique.} 

The divergence of (T>) at nonzero T indicates the breakdown of our approximation. 
We may see this more explicitly by looking at the argument of the exponential in (J35j) . 
The Taylor series J2^=o jt for e x is dominated by the terms for which j is O(x). At T = 0, 
the argument of the exponential appearing in (J3o"|) is . f?w dpj 6 = fog^^), which 
is (9(logL). Therefore, the dominant contributions to the exponent come from values 
of j that are likewise 0(\ogL). However, when T > 0, the average number of quanta 
present in each of the lowest-lying modes is 0{V) unless (3 > (fFPmin) -1 - The argument 
of the exponential and the dominant values of j are then O(LlogL). This contradicts 
our assumption that j is smaller than 0{L). 

We have shown that when the temperature is nonzero, we must receive a substantial 
contribution from terms for which j is O(L) or larger, so the methods we have used 
so far do not apply to this case. However, the T = treatment does provide us with 
several insights that will prove useful when we analyze the case in which T > but the 
"transverse temperature" is still vanishing. We shall return to that case later; however, 
there are several more points about the zero-temperature situation that we must consider 
first. 

Thus far, we have assumed that a -1 is the correct cutoff for the momentum integral 
appearing in ()35|). However, there is also another cutoff we must consider. In order to 
obtain a physically meaningful result, we must determine which cutoff is the relevant one. 

A real Luttinger liquid has a finite bandwidth vpa -1 , because the phonon description 
breaks down at large values of Pi. The Luttinger model only describes the low-energy exci- 
tations of the fermion system. Whenp^ becomes comparable to hp, the energy — 
of an excitation can no longer be approximated by vppi and is no longer much smaller 
than ep, so we expect the physical cutoff a -1 to be 0{kp). 

The second cutoff arises from the energy denominator (J31)) . When pi is large compared 
to vmA, the integrand decays as p~ 3 . If there existed an intermediate range of pi values, 
where mA/kp <C Pi <C VmA, the integrand would decay as p~ 2 in this region. However, 
since we must have A > -^-kp in order to have an energy gap at all, such a region can 
not exist. So the rapid p~ decay at large values of pi will provide a cutoff for the integral 
which we expect to be 0(V mA). 

We must determine which of these two cutoff scales is smaller, since the smaller one 
determines the behavior of the energy shift. When A ^> ep, the first cutoff, near the 
Fermi momentum, is smaller. However, when the energy gap is comparable to the Fermi 
energy, \pmA is itself Oikp). Without more precise values for the cutoffs, it is not clear 
which one matters in this situation. 

We shall now estimate the cutoff p' max arising from the p~ 3 behavior of (jSH). To do 
this, we shall again make use of the fact that the integral over the pi is dominated by 
the region where the pi are small. The value of p' max is determined by the behavior of 
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the integrand when not all the Pi are small; at least one of the pi must be comparable 
to VmA. However, the dominant contribution under such conditions comes from the 
region where only a single Pi is large compared to p mm - If more than one Pi is large, the 
contribution is much smaller. This means that we may determine p' max by analyzing the 
case in which all but one of the Pi are negligible. 

More precisely, we approximate b k ^ F+ Yf_ Tn b y S k ,k F - Pl and £i=i ±u Pi by u Pi . The 
choice of signs corresponds to the fact that the virtual phonon at momentum pi must be 
created, not annihilated, since almost all the phonon modes are empty of quanta when 
T = 0. One of the Pi integrations is now nontrivial. We shall perform this integration 
explicitly and extract the cutoff from it. Our prescription for doing this is that 



dpi 1 



Pi a + ^[(k F - Pt y - e F 

The integral in (|4H|) is elementary. The result is 

50 dpi 1 




(43) 



2m 




2mA + e^p 2 min , (44) 



so, neglecting e 27 p min compared to 2mA, we have that p' max = e 7 v2mA. (This is really 
the appropriate value of p' max only when q <C v2mA. A more careful calculation, which 
takes into account the g-dependence of the cutoff, is located in Appendix iBl) 

Although this estimate is rather inexact, it shows the general character of the cutoff. 
At A = ep, the cutoff takes on its minimum value of p max ~ 1.78k F . However, we have 
not yet determined whether this is ever the relevant cutoff. 

We shall defer the final resolution of this question until Section 0J In the meantime, 
we shall make the replacement 

+00 p -api / r „ x 

dpi ► log ( = ) , (45) 



el 



Pi V 27T 



where in all cases, p max is 0{kp). Using the prescription (|43j) to evaluate the integral and 
inserting our expression for V(T = 0), we get 



1 



The denominator, A + -^Aikp — sq) — k F ] + su g = A + ^— is independent of s, so we may 
trivially perform the sum over s. Thus far, we have restricted our attention to the case 
for which c k corresponds to a right-moving fermion. The terms in which c k corresponds 
to a left-moving particle also contributes an equal amount, so we have 

.„ |gM| 2 <HW / 1\ , 47 , 
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(14 7j) is the energy shift for a single bosonic mode, with momentum +q. To obtain the 
total shift for the system, we must add in the similar renormalization of the left-moving 
mode with momentum — q and integrate over q. Since the density of momentum states is 
proportional to L, the total shift is extensive, as it must be. We must also sum over all 
the possible n p values to which dk may correspond; this will give an infinite sum of terms 
of the form (J4~T|) . in which A takes on all the values Ao,n p , n p > 0. 

We shall now turn our attention to the T > case. The crucial simplification that 
allowed us to evaluate the T = energy shift in closed form was the fact that the energy 
defects were effectively independent of the intermediate states. A remarkable cancellation 
then occurred between the prefactor V and the exponential that arose when we summed 
over all all intermediate states. A similar cancellation will occur here. 

When T > 0, we may not necessarily neglect Y?%=\ TPi and YLl=\ ±u Pi , because j is 
not generally small. The presence of these two sums gives rise to a nontrivial energy 
denominator. However, if A is formally very large, so that any other energy may be 
neglected compared to it, the energy defect again becomes independent of the Pi. As 
in the previous case, the temperature is negligible when compared with the transverse 
energy scale; that is, the "transverse temperature" is effectively vanishing. With this 
simplification, it will again be possible to obtain an analytical expression for the energy 
shift. 

If we again assemble the various terms that compose the energy shift (now assuming 
that A is very large), we have 



AE.L 



E 

s=± 



i \g(q)\ 2 q(n+ q + 



2 



x 27T« 



2tt 



E 

{m +p } 



({m +P }\ n p >oe A+p <- A +^+1W P }> 



A 



(48) 

where X +p = —^2%/ (pL)e~ ap ^ 2 , and A\ p and A +p are the canonical ladder operators 
corresponding to the p + (±p) quanta. 

Because the energy defect in (|48|l is independent of the {m +p }, we may rewrite the 
energy shift as 



AE, 



\g(q)\ 2 q(n +q + I / *_ x 

E Tfa-T ~^ n ^\ II [e +P +P +P +P 

s=±l ^ n I p>0 v 



X 



E \{ m + P }}({ m +p}\ 

{m +p } 



n 

p>0 



e \ +p A< +p -\* +p A +p 



IKp})- 



(49) 



The operator in brackets is just the identity, so we may drop it. Furthermore, since 

t r / m-1 



A 



+P A ] +p -\* +p A +p is antihermitian, n p >o ( e x+pA +p~ x *+p A+p 



n 



p>0 



-A* A A 



and the entire matrix element expression is unity. [We now see that the remarkable sim- 
plicity of the T = result is a consequence of the fact that when the energy denominator 
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is constant, the sum over intermediate states may be evaluated simply using the closure 
relation. However, the detailed analysis that led up to was necessary for us to justify 
taking the energy denominator to be constant.] 

This T > result shows that, when A is very large, we recover a similar expression 
to the one we obtained in a different extreme limit, that of zero temperature. The total 
energy shift is 



This expression is actually identical to the large A limit of (I47|) . In that limit, may be 
neglected compared to A, and the correct cutoff p max is just a -1 , because p' max grows as 
O(ymA). However, we still do not know to what actual value a corresponds. We shall 
now turn our attention to a different model, which will allow us to answer that question. 

4 Energy shift for a more physical model 

In this section, we present an alternative model — one which is useful when T = and 
A is very large. Using this model, we shall again calculate the energy shift. From our 
calculation we may extract the physically relevant value of a -1 , the cutoff arising from 
the breakdown of the phonon description. Moreover, we shall find in Section El that our 
calculation requires only slight modification in the presence of strong interactions between 
the low-lying fermion states. 

The Luttinger model does not accurately represent the fermions lying deep within the 
Fermi sea, so it is not necessary to use a bosonized form for the matrix elements of Ck 
when k 96 hp. Since a basic assumption of the Luttinger liquid theory is that the presence 
or absence of additional fermions far from the Fermi surface should have minimal effect 
on the low-energy excitations, we shall simply assume that the action of Ck (k ^ kp) on 
a Luttinger liquid state is to create a hole in the Fermi sea, without affecting the phonon 
configuration in any way. The only exception to this prescription is that when the fermion 
operators are combined in the form (0), as they are in H 2 , we shall used the bosonized 
expressions. 

Our approximation will require that A ^> ep. This condition ensures that the energy 
shift is not dominated by the terms with the smallest energy defects; instead, the defects 
are dominated by A and depend only weakly on the energy of the fermion annihilated by 
Cfc. So the entire Fermi sea will contribute to the energy shift, and the region near the 
Fermi surface (where our approximation for Ck is invalid) will make a comparatively small 
contribution. The condition that T = ensures that the occupied states are precisely 
those with \k\ < kp. 

Now we must identify the states between which we shall be calculating matrix elements. 
As before, we shall have the state |{^+ p }) on the right. We shall find the energy shift of 
this state due to interactions with the states \{n' + }; k^q;k). These states have bosonic 




(50) 



,12 
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occupation numbers {n' + } which are the same as {n +p }, except that n'_ 



+q 



+ 1. 



These states also have an excited fermion with momentum k =p q and a hole in the Fermi 
sea with momentum —k. 

The matrix elements we shall need are then 



(Kb,}; k - sq] k\p+(sq)d\_ sq c k \{n +p \) = (n +q + ^-)0(k F - \k\) 



2 J 



(51) 



for s ± 1. The step function ensures that the matrix elements vanish when \k\ > kp. The 
energy defects are 



E{n'};k-sq;k ~ E { n+p } 



2m 



(k — sq) 2 — k 2 + sujq, 



(52) 



so the energy shift is 
AE = 



— I ( )\ 2 V V KK+p)' k ~ s ^ k \P+( S( l) d k-s g c k\{n+p})\ : 
L 2 Wm k h F ski A + ^ [(k - sqf - P] + su, g 



n 



+q 



+ 



5 + 1 



1 A + ±[(k-s q y-k*]+su q 



(53) 



Converting the sum over k into an integral, we have 

8 + 1 



AE 




—) r dk^- 

2nJ J-k F (A + 



2m 



+ SUJ a ) — s — k 



A + £- + su a )-s±k 



2m 

Isqkpjm 
A~+f- 

2m 



(54) 



We have already required that A be large compared with the Fermi energy. Since 
q <C kp, we may use the approximation log(l + x) ~ x, approximating the logarithm in 
(|54|) by its leading term. (In fact, because of cancellation between the s = ±1 terms, 
this is equivalent to approximating the logarithm by its two leading terms for purposes 
of determining the frequency.) This gives us 



AE 



l#(<?)| 2 #i 



7T Z 



A + f- 

2m 



n. 



~ q+ 2 



(55) 



As in Section El we must sum over all values of q to get an extensive energy shift. 

We now compare the result (J53)l with (|47j). The models leading to these expressions 
are both valid when T = and the energy gap is large, so the two expressions must agree 
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under these circumstances. In fact, ()47|) and (JS3j) agree in their region of common validity 
exactly if p max = k F . Since k F is substantially smaller than any possible value of p' max (for 
large A), we must conclude that the relevant cutoff is p max = a -1 . The value of a~ l is a 
property solely of the Luttinger liquid, and a~ l is thus idependent of A. Therefore, a -1 
is always equal the Fermi momentum when T = 0. We then see that a~ l < p' max for all 
values of A, so a~ l is always the most relevant cutoff, and (|B3j) is the general expression 
for the T = energy shift. (We shall discuss the possible relationship between a and the 
temperature in Sectional) 

A simple estimate of the magnitude of any correction to (|53|) due to p' max is that the 
fractional error should be approximately 0{k F /p' max )- Since we found a minimum value 
of Pmax °f only 1-78/cf, this error could be significant when the energy gap is small. This 
must be kept in mind when this formula is applied to systems for which A is 0(€f)- 



5 Extensions and discussion 



The model outlined in Section 0] may be applied to systems with more general interac- 
tions. In particular, our techniques may be applied to the interacting Tomonaga-Luttinger 
model. The Tomonaga-Luttinger Hamiltonian Ho + H\ is diagonalized by the Bogoliubov- 
transformed operators [21] 



(56) 



where 



X 



VF + sM + sM 

Z7T Z7T 



1/4 



n„ _i_ 34(g) 



32(g) 

2tt 



(57) 



The p± satisfy the same commutation relations as the p±. To determine the effects of the 
additional interaction H 2 , we only need relate p+(q) + p~(q) to p'Aq) + p'-{q)- From ()5b|). 
it is obvious that 



p+(q) +p-(q) 



32(g) 

2w 



1/4 



Z7T Z7T 



[p'M + p'M)}, 



(58) 



so the calculation of AE requires only two modifications. We must use the Tomonaga- 
Luttinger frequency 



to „ 



9 



94 

Ztt 



92_ 

2tt 



n l/2 



(59) 



in place of u q , and we must multiply the whole expression for the energy shift by |X| 2 . 
This complicates the expressions somewhat; however, if we neglect terms of 0(u q ) com- 
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pared to A (since q <C kp), we simply have 

1/2 i 

+ 2) • ( 6 °) 

It is slightly more difficult to perform this calculation using the bosonized operators. The 
mixing of the the right-moving and left-moving phonons in ip± adds to the complexity of 
the calculations El 123 I2H] ■ Performing the calculation, we find that the matrix ele- 
ments factorize into p' + terms and p'_ terms. The cancellations that appeared in Section [3] 
at T = now occur separately for the two types of terms. The end result is the same as 
flHOjl . provided we set p max = or 1 = kp and again neglect the 0(w q ) terms. 

Our main result is l)55j) . the expressions for the energy shift. It is interesting to note 
that remains regular — without a pole or a branch point — at A = e F . When A 01 = e F , 
we are in a slightly different regime from the one we have considered, because some n p = 0, 
£ = ±1 states of the trap will be occupied. However, we can still draw some conclusions 
from the behavior of the energy shift. Although the ground state changes qualitatively 
at this point, as fermions "spill over" out of the Luttinger liquid into the radially excited 
state, the phonon energies remain finite. Even if the decay of a phonon into a n p = 0, 
£ = hole and a n p — 1, I — fermion is allowed energetically, the phonons appear to 
remain metastable against this decay channel. This represents an interesting "softening" 
of the expected phase transition at this point. However, in order to understand this 
situation properly, other interactions than (jl(J|) . including those involving the occupied 
£ = ±1 states, must also be considered. 

Since the cutoff parameter a is related to the renormalization group flow for the 
Luttinger liquid [21], the identification of kp with a~ x has interesting implications. In the 
Luttinger model, a is effectively a free parameter; its value needs to be specified on the 
basis of physical considerations external to the model. Our results at T = suggest that 
vpkp = 2ep is the maximum physically meaningful bandwidth for the system. This order 
of magnitude for the cutoff has a clear physical basis, as discussed in Section [3] Since 
the Luttinger model (with spin) is known to be equivalent to several other condensed 
matter models [2D 123 I2E1 123 EDI EH and field theory models H3E21, the specification of 
a particular value of a may have interesting implications for the related cutoffs of these 
analogous systems. 

Although we have identified a -1 = kp as the correct physical cutoff at T = 0, our 
results do not tell us anything about the renormalization group flow for nonzero temper- 
atures. Moreover, although (joT)|) is valid for finite temperatures, it is only meaningful in 
an extremely singluar limit; all the temperature-dependence of the expression has been 
removed by making A very large. The temperature-independent expression (|5U|) is consis- 
tent with a nontrivial T-dependence of the energy shift, becuase AE is actually formally 
vanishing in this limit, since it is proportional to the ratio [|5 , (<?)| 2 ct~ 1 <?] /A. So the ques- 
tions of the finite-temperature energy shift and the nature of the renormalization group 
flow for T > remain open. 



AE 



7T 2 A 



„. I 94(g) 92{q) 

V F + ~2i 2lT 



Vp + sM + 

Z7T Z7T 
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We have used several complementary techniques to evaluate the energy shift. The 
bosonized method used in Section |3] is valid in a wider range of situations; however, we 
needed the calculations of Section HI to determine the unknown constant p max appearing 
the bosonized result. We found a frequency shift that is always 0[\g(q)\ 2 qkp/ A], even for 
nonzero temperatures or in the presence of strong Tomonaga-Luttinger-type interactions. 
When experimenters are able to build fermionic traps with sufficiently high aspect ratios 
to observe Luttinger liquid behavior and to trap fermions at sufficiently low temperatures, 
formulas such as these should become readily testable. 
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A Appendix: Harmonic oscillator matrix elements 

To obtain the general formula (|20|1 . we first apply the Baker- Campbell-Hausdorff formula 
to the operator e XAt ~ x * A . Since [— \*A, \A^] = — |A| 2 commutes with both A and A\ we 

e AAt-A*A = e -A* V At e |A|72. (61) 

This reduces the problem to the determination of (n\e~ x * A e XA \m). 

There are three cases: m = n, m > n, and m < n. The first case is the simplest. 
If m = n, each factor of A from e~ x * A must be paired with a factor of A^ from e Aj4t ; 
otherwise, the term does not contribute. So we get 



(m\e- x * A e XAt \m 



) = (m|l + (-AM)(AAt) + ^(-AM) 2 (AAt) 2 + ...|m) 



fX-D'lAi* 1 (m + 2)! 



i=0 



(i\) 2 m\ 



= F(m+l;l;-|A| 2 ) . (62) 

For the case of m > n, there must be I more factors of A^ than factors of A for a term 
to contribute. This immediately leads to the series 

He^VV) = (m|i(A^)' + ^(-AM)(A^) ,+1 

+ TO ( - AM)2 ( A < V " |n) 



2i [ m + % ) 



l\ toi\(l + i)r 711 ml 
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4v^ F ( m+i;,+i; - |A|2 )' 



(63) 



In the third case, I is negative. The calculation proceeds along essentially the same 
lines as for the I > case. There must be |/| more factors of A than A* for a term to be 
nonzero, so we get 

{m\e 



\*a 6 \a^ 



In) = (m\-^(-X*A)^ + ^jyp— q? (-X*A)^ +l (XA^ 
^^(-A^) W+2 (AAt) 2 + -|n) 



;-i)»|A 



i,M2i(n + «) ! H 



-A*) 



1 n\ 




F(n + l;|Z|+l;-|A| 



n\ v m\ 

2 s 



(64) 



Combining equations f2J), (JH2J), (p^jl . and (fBlj) . we get CT . 



B Appendix: Momentum dependence of p f 



max 



We shall now examine how the cutoff p' max depends upon q. We do this by modifying the 
prescription (}4*3j) for p' max to 



f°° i « 1 w ^= (65) 

To do this, we use the integration formula 



+oo 



£nsp api + fyo 2 + cpf a \ p min V c 




a + be Tp min + ce 2 Tp^ 



(66) 



b . / 2ic \ / — 2ic 

i loe 



2a\J-b 2 + 4ac [ \V~-b 2 + 4ac / + 4ac/ 

/ V — ^ 2 + 4ac — ib — 2ice~ 1 p min \ ( \J — b 2 + Aac + ib + 2ice~ 7 p I1 

+2 log ===== -i log 



+ 4ac / V \/-b 2 + 4ac 

In this case, we have a = A + b = — , and c = We shall define x to be the 

2m ' m ' 2m 

frequently- appearing ratio x = — ^_ 62+4ac = ~ s ^La " ^ e sna ^ ^ nc i that Pmax depends 
upon 5 only through x. 

To simplify the expression ()66j) . we first note that since p m i n is small compared to 
any other momentum scale, we may approximate a + be~' y p min + ce^p 2 ^ ~ a and ib + 
2ice~ 1 p m i Il ~ ib. This reduces the first term on the right-hand side to - log {-^— \ -)■ 
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We now turn our attention to the complex terms. Since these terms appear in conju- 
gate pairs, we need only calculate their real parts. Each term has an overall factor of i 
times a logarithm, so we shall only need the logarithms' imaginary parts. Because ^_^ 4ae 
is purely imaginary, we may simplify the second and third terms on the right-hand side 
of (JUD to 



x 
2a 



i log 



2ic 



\J-b 2 + Aac, 



— i lot 



-2ic 



V-b 2 + Aac / 



x 



~2a ^ l0g ^ ~ ilo s( -i ) 



x 

~2a 

TlX 

2a~' 



71 



71 



(67) 



We may evaluate the last two terms of (|B1)|) similarly, getting 



x 



[log(l + ix) — i log(l — ix)] 



So the entire expression becomes 
'"+ 00 dpi 



x . 

~2a 

— tan -1 x. 
a 



ztan l x — itan 1 (—x) 



api + bp 2 + cpf 



1 , fe^Ly/2mA + q' 2 \ 1 
log — 



2tt 



( -i 1X \ 
x tan x H — 

. V 2) 



(68) 



(69) 



We may absorb the second term on the right-hand side of ()69|) into the logarithm by 
exponentiating it. We then identify the cutoff as 



Pmax = e 7 -\/ 2mA Vl + x 2 exp 



/ -1 ^ 
x tan x H — 

V 2 



(70) 



This expression depends upon s. However, only the value of p' max averaged over s = ±1 
is ultimately relevant, so we may modify our expression for p' m&x to become 



Pmax = e 7 v / 2mAv / l + x 2 exp ^xtan 1 x^j cosh ^"^"j • (71) 
For fixed A, this takes on its minimum value in the q = case considered previously. 
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